How Flat is Our Universe Really? 
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Distance measurement provide no constraints on curvature independent of assumptions about the 
dark energy, raising the question, how flat is our Universe if we make no such assumptions? Allowing 
for general evolution of the dark energy equation of state with 20 free parameters that are allowed 
to cross the phantom divide, w(z) = 1, we show that while it is indeed possible to match the first 
peak in the Cosmic Microwave Background with non-flat models and arbitrary Hubble constant, 
H , the full WMAP7 and supernova data alone imply -0.12 < Q k < 0.01(2cr). If we add the HST 
Ho prior, this tightens significantly to = 0.002 ± 0.009. These constitute the most conservative 
and model-independent constraints on curvature available today, and illustrate that the curvature- 
dynamics degeneracy is broken by current data, with a key role played by the Integrated Sachs 
Wolfe effect rather than the distance to the surface of last scattering. If one imposes a quintessence 
prior on the dark energy (— 1 < w(z) < 1) then just the WMAP7 and supernova data alone force 
the Universe to near flatness: fl^ — 0.013 ± 0.012. Finally, allowing for curvature, we find that all 
datasets are consistent with a Harrison-Zel'dovich spectral index, n s = 1, at 2a. 



Introduction - The sign and magnitude of the cosmic 
curvature, ilk, is one of the most fundamental charac- 
teristics of our cosmos. The sign controls the topology 
of the universe while the magnitude has real importance 
in testing theories: eternal inflation would be seriously 
tested if \ilk\ > 10~ 4 [I, 2] while anthropic considera- 
tions that suggest that \Slk\ might not be vanishingly 
small [3] . Assuming ACDM we are not far from reaching 
this regime with the latest curvature constraints around 
the <7fi fc ~ 10 -3 level [4, 5]; a huge improvement over the 
major break-through from the BOOMERANG mission a 
decade ago which gave \ilk\ < 0.2 [6], itself an order of 
magnitude improvement over earlier constraints [7]. 

The key step forward in measuring curvature was the 
realisation that the position of the first acoustic peak 
in the Cosmic Microwave Background (CMB) provides a 
standard ruler, and hence the distance to the surface of 
last scattering [8]. However, this alone does not constrain 
flk because of the well-known geometric degeneracy be- 
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tween ilk and the Hubble parameter today, Hq [9] . With 
the addition of current large scale structure data, how- 
ever, this degeneracy is now almost completely broken. 

We are concerned here, however, with the more perni- 
cious degeneracy between ilk and dark energy dynamics, 
parametrised through the equation of state, w(z), which 
even an infinite number of perfect distance measurements 
cannot break [10-15]. It follows from simple degrees of 
freedom counting - we are trying to measure a free func- 
tion - the expansion rate, H(z) - as well as a constant, 
ilk, from a single free function, the distance dA(z). Such 
degeneracies are not peculiar to distances alone [15], nor 
to curvature, as illustrated by the dark matter-dark en- 
ergy degeneracy which makes it impossible to measure 
f2 TO without assumptions about dark energy [16]. 

The dark energy equation of state, w{z), which will 
make a curved universe with curvature parameter ilk ex- 
actly mimic the distances in a flat, ACDM model at all 
rcdshifts, is given by [13]: 
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where Dl = (H$/c)dL is the dimensionlcss luminosity 
distance in the flat ACDM model and primes denote red- 
shift derivatives. Fig. (2) gives the CMB power spectrum 
for an open model with = 0.15 with w(z) given by Eq. 
(1) showing that the it matches the first peak, and hence 
the distance to the surface of last scattering, perfectly. 



Only by making assumptions about the dark energy 
dynamics can any distance measurements, even if mea- 
sured precisely at all redshifts, provide constraints on 
curvature. However, such constraints are fictitious un- 
less we know the true nature of dark energy, which is not 
the case today. The strongest constraints arise, of course, 
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FIG. 1: 100 randomly-selected dark energy w(z) curves from 
the converged chains with Prior I ( — 15 < Wi < 1) and the 
BASIC data package (WMAP7+SNIa+BBN). Darker curves 
have higher likelihood. While w(z) is well-constrained at very 
low redshift, it is unconstrained at z > 1, filling the entire 
prior. 



by assuming that the dark energy is a cosmological con- 
stant, A with w = — 1. 
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FIG. 2: Temperature angular power spectra for various 
curved models with the binned WMAP7 data. Blue solid 
curve: Ct for the best fitting model with VLu = —0.03 and 
H = 63.4 kms~ 1 Mpc~ 1 . Red thick dash-dotted curve: The 
best-fit from an MCMC search with fixed flk = 0.15. It 
has Ho = 56.4 kms _1 Mpc _1 . The fit is essentially perfect 
for £ > 10 but the overall fit is poor due to the large In- 
tegrated Sachs Wolfe effect which disfavors all significantly 
open models. Red dashed curve: With fixed fife = 0.15 
and Ho = 71 kms _1 Mpc _1 , we use w(z) from Eq. (1) that 
matches the flat ACDM distances at all redshifts. Although 
the first peak matches perfectly, the remaining fit is very poor. 
These two examples show that the distance to the surface of 
last scattering is not the primary driver of our constraints on 
curvature. 



Since essentially all current constraints on curvature 
derive from distance measurements from the CMB, su- 
pernovae and Baryon Acoustic Oscillations (BAO) with 
restrictive assumptions about the dark energy, a natural 
question is 'how curved is our Universe allowed to be if 
we do not assume anything about the dark energy dy- 
namics'? Could we fit the data with 51^ = 1.5 with a 
suitable w(z) for example? 

There are two relevant issues to answering this. The 
first is the freedom in the w(z) parametrisation: how 
many free Wi parameters? The second is the range over 
which the Wi parameters are allowed to vary. To match 
fiat ACDM distances with a closed model requires aw(z) 
that typically must go out of the range — 1 < w{z) < 1 
allowed by the weak energy condition and naive non- 
superluminal speed of sound [13]. The w(z) freedom dic- 
tates how subdominant dark energy was to matter and 
radiation when these components dominate the universe. 
Without enough w(z) freedom a distance measurement 
during matter domination in addition to the CMB docs 
break the curvature degeneracy [17]. 

Method - We used the CosmoMC Markov Chain Monte 
Carlo (MCMC) package [18] together with the WMAP7 
likelihood for our parameter estimation and typically 
used five chains of about 3 x 10 5 steps each and used 
the Gelman-Rubin test for convergence. To allow for 
general dark energy evolution, we used the PPF module 
[19] which splines an arbitrary number of parameters, 
Wi, at arbitrary redshifts to give w{z). We found that 
our results were not very sensitive to the total number 
of Wi parameters used which is not surprising since the 
effective number of accurate distance measurements is 
relatively small. We simply need to ensure that we have 
enough degrees of freedom for w(z) to capture the de- 
generacy. All results in this paper are quoted for 20 Wi 
spline coefficients located at logarithmically-spaced red- 
shifts between z = and z = 6000, allowing for rapid 
w(z) variation at low redshift. 

We considered two classes of dark energy priors on 
the Wi. Prior I was flat over the range iu, € [—15,1] 
[37] while Prior II was flat over [—1,1], mimicking 
quintessence-type models [13]. We assume adiabatic 
initial conditions and used fco = 0.002 Mpc -1 as the 
pivot scale for scalar and tensor power spectra nor- 
malization [29]. The parameter set we used in our 
MCMC chains, which includes the 9 standard cosmo- 
logical parameters in addition to our 20 Wi, is: V = 
(n b h 2 ,n c h 2 , [9 or H ],T,n k ,n s ,A s ,r,A sz ,Wi). 

We used two minimal sets of data that are chosen to 
be maximally independent of assumptions about curva- 
ture and dark energy. The BASIC data package con- 
sists of the CMB temperature and polarization data from 
WMAP7 [20], the Union2 supernovae [21] and the Big 
Bang Nucleosynthesis (BBN) prior VL b h 2 = 0.020 ±0.002 
(2a) [22]. In addition we consider the BASIC + HST 
data package, which adds the Hubble Space Telescope 
constraint [23] on the Hubble constant, H = 72 ± 
8 kms _1 Mpc _1 [38]. We do not include BAO and growth 
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FIG. 3: Joint contours on ilk and n 3 for the three data sets 
and priors: BASIC + Prior I -15 < Wi < 1, BASIC + 
HST + Prior I, BASIC + Prior II (-1 < Wi < 1). Here 
BASIC = WMAP7 + SNIa + BBN. Even allowing for general 
dark energy models and conservative data choice constrains 
the curvature to satisfy ilk = 0.002 ± 0.009, but all cases are 
compatible with n 3 = 1. 



information to be conservative, since the assumption of 
flatness and ACDM comes in subtly, for example in the 
N-body simulations used to calibrate results. As we now 
show, we do not need extra data to rule out significant 
curvature even in the case of the most general w(z). 

Results - In Fig. (1) we show a random selection 
of 100 splined w(z) curves from our chains with BA- 
SIC data with higher likelihood curves shown, highlight- 
ing the freedom accorded to the dark energy in our anal- 
ysis. Three relevant CMB spectra arc shown in Fig. (2). 
The global best-fit with Prior I and the BASIC data 
package which has ilk — —0.03 is shown as the solid 
blue curve. The best fit with fixed ilk = 0.15 and 
i?o = 71fc kms~ 1 Mpc _1 (red thick dash-dot curve) shows 
why open models are ruled out by the Integrated Sachs 
Wolfe (ISW) effect at I < 20. This is not surprising 
since to match flat ACDM distances typically requires 
very rapid w(z) evolution at low redshifts [13] with re- 
sulting large ISW effect [26-28]. The C e 's of the model 
with fixed ilk = 0.15, WMAP7 best-fit parameters and 
the w(z) corresponding to Eq. (1) which ensures that 
the distances are identical to the best-fitting flat ACDM 
model at all redshifts is shown by the dashed line. Clearly 
the distance to the surface of last scattering plays little 
direct role in constraining curvature in general, unlike 
the case for ACDM. 

Our main finding is that with the BASIC + HST data 
we recover a constraint on the curvature: ilk = 0.002 ± 
0.009, even with effectively no limits on the Wi (Prior 
I). The HST prior on the Hubble constant is critical in 
removing the closed branch of universes that are excellent 



fits to the BASIC data package and have low Hq. Indeed, 
without the HST prior the ilk posterior peaks around 
ilk = —0.085. Assuming quintessence-like dark energy 
(Prior II) also removes the closed branch, irrespective of 
the HST constraint. The importance of allowing crossing 
of the phantom divide, w{z) = —1, is a consequence of 
trying to match distances, since Eq. (1) requires w — >■ 
— oo at some redshift when ilk < 0. These results are 
exemplified in Fig. (3) and summarised in Table (I). 

It is a subtle combination of effects, wide coverage and 
multiple datasets, that gives good constraints on curva- 
ture in the case of general dark energy dynamics. One 
of the main reasons for this result is that while it is pos- 
sible to match either the flat ACDM distances or ex- 
pansion rate, H{z), in a curved cosmos, it is not pos- 
sible to do both simultaneously over an extended red- 
shift range: ilk = is the only solution to the equation 
sm(y/-flkx)/\ / -Qk = X, where \ = H J dz'/H(z') is 
the usual flat ACDM comoving distance. Hence, simul- 
taneous low-redshift measurements of both Type la su- 
pernovae (SNIa) which constrain distances and the ISW 
effect, which constrains H{z), provide constraints on cur- 
vature even in the absence of an H prior. 

In addition to our BASIC and BASIC + HST runs, 
we also undertook runs with the CMB lensing results 
from the Atacama Cosmology Telescope (ACT) measure- 
ments [31] and the time-delay distance to the lens system 
B1608+656 at z = 0.63 [32] as alternatives to the HST 
prior. However these two additional datasets did not sig- 
nificantly impact the constraint on ilk- The impact of 
ACT is shown in the final two columns of Table (I). 

We can also look at the effect of general dark energy 
dynamics and curvature on other parameters, such as the 
scalar spectral index, n s , where assuming flat ACDM, 
there is evidence for n s < 1 at a significance of about 
3a [33]. As shown in Fig. (3) our results are all consis- 
tent with the Harrison- Zel'dovich value of n s = 1, due 
to the larger ISW effects from both curvature and dark 
energy dynamics [28], illustrating how cosmological con- 
straints on the early universe are tightly coupled to the 
late universe. 

Conclusions - We have established constraints on the 
curvature allowing for a general evolution of the dark 
energy equation of state which does not rely on ar- 
tificial breaking of the curvature-dynamics degeneracy. 
The WMAP7 and Union2 supernova data are sufficient 
to constrain the curvature to ilk = 0.013 ± 0.012 if 
quintessence-like behaviour are assumed but if crossing 
of the phantom divide w(z) = — 1 is allowed then we 
find -0.12 < il k < 0.01 at 2cr. With the HST prior on 
H added we find il k = 0.002 ± 0.009, implying that the 
curvature-dynamics degeneracy does not change the re- 
sults of the standard ACDM picture although we empha- 
sise that our curvature constraints are no longer primarily 
driven by the distance to the surface of last scattering, 
as illustrated in Fig. (2), but rather the full shape of the 
power spectrum at both small and large angular scales. 

While these results do not assume that we know the dy- 
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0.63 ± 0.06 


0.75 ± 0.02 
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TABLE I: Summary of results. WMAP7 refers to the results in [29] where only the WMAP 7-year data are used. OACDM has 
f2fc ^ 0. All other runs use 20 Wi parameters with Prior I (II) imposing —15 < Wi < 1 (—1 < Wi < 1). The BASIC data consists 
of WMAP7+SNIa+BBN while HST is the Hubble constant constraint H = 72 ± 8fcms~ 1 Mpc- 1 and ACT is the CMB weak lensing 
spectrum measurement from ACT [30]. Note the large range of ages in cases where the HST prior is not used and that n s = 1 is compatible 
at 2a with all dynamical dark energy runs. 



namics of the dark energy, they do, however, assume that 
gravity is described by General Relativity and the speed 
of sound of the dark energy is unity. In this sense our 
results still leave a last loop-hole. Fortunately number 
counts, N(z), time delays, lensing and growth of struc- 
ture measurements depend on curvature very differently 
compared with distances and hence future data will allow 
us to close this loop-hole too [11, 15, 34, 35]. The key 
will be to extract accurate measurements of these observ- 



ables that are truly independent both of the assumptions 
of flatness, ACDM, and General Relativity. 
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